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There is very limited knowledge of the kinematical relations for the velocity structure functions 
higher than three. Instead, the dynamical equations for the structure functions of the velocity incre- 
ment are obtained from the Navier Stokes equation under the assumption of the local homogeneity 
and isotropy. These equations contain the correlation between the velocity and pressure gradient 
increments, which is very difficult to know theoretically and experimentally. We have examined 
these dynamical relations by using direct numerical simulation data at very high resolution at large 
Reynolds numbers, and found that the contribution of the pressure term is important to the dy- 
namics of the longitudinal velocity with large amplitudes. The pressure term is examined from the 
view point of the conditional average and the role of the pressure term in the turbulence dynamics 
is discussed. 

PACS numbers: 47.27.Ak, 47.27.Jv, 47.27.Gs, 05.20.Jj 
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INTRODUCTION 

In Kolmogorov's theory (1941, hereafter K41), the mo- 
ments of velocity increments in the isotropic turbulence 
scale as (U p ) cx r^ p with ( p = p/3, where U is the lon- 
gitudinal increment defined as U = (u(x + r) — u(x)) ■ 
r/r[l, 2]. However, it is now widely accepted from the ex- 
perimental and direct numerical simulation (DNS) data 
that at high Reynolds numbers the intcrmittency builds 
up with decrease of scales and scaling exponent C, p does 
not obey K41 and is a non-decreasing function of p. 

Arguments have recently been raised as to whether or 
not the scaling exponent of the longitudinal structure 
function is equal to that of the transverse one at the 
order higher than three [3-9]. The structure function 
(moment) of the velocity increment is defined by 



Stokes equation as 



S m , n (r) = (U m V n ), 



(1) 



where V is the transverse velocity increment. In the in- 
ertial range at very high Reynolds numbers, S m ^ n (r) is 
assumed to obey power law S m . n (r) oc r^ m,n . In the ho- 
mogeneous isotropic turbulence, the incompressible con- 
dition yields the well known relations at the second and 
third order 



SoAr) = S 2 ,o(r) + T^ 5 " 2 ' ' 

Si, 2 (r) = is 3 ,o(r) + ~S 3 ,o(r), 
6 6 dr 



from which it follows £2,1 



Co, 2 and £ 



3.0 



Ci 



(2) 

(3) 
Kol- 



S 3,o( r ) 



S 3 ,o(r) + 25i, 2 (r) 



4_ 



(4) 



so that £3.0 = C1.2 = 1, where e is the average rate of 
energy dissipation rate per unit mass. However, for the 
order higher than three, no equations are known that 
tie directly the longitudinal structure functions to the 
transverse or mixed structure functions. Instead, rela- 
tions among various type of structure functions are de- 
rived from the Navier-Stokes equation which contains the 
pressure and viscous terms [14, 18, 20] such as 



<(<5 r (V||p))W), ({8 r {V 2 U))UVV«) 



where 



dp(xi) dp(x 2 ) 



11 \ 0x1 0x2 
S r (V 2 U) = (V 2 Xl u( Xl ) - W 2 X2 u(x 2 )) 



(5) 

(6) 
(7) 



mogorov's 4/5 and 4/3 laws are derived from the Navier- 



Gotoh considered the pressure term from the view 
point of the conditional average of the longitudinal in- 
crements of the pressure gradients (5V|i£>|t/, r) for given 
value of U (see later sections for more detail) [21]. The 
conditional average is expressed in terms of the quadratic 
function of U and the coefficients are fixed by some sta- 
tistical constraints. The qualitative feature of the con- 
ditional average was compared with DNS data and the 
results was found to be satisfactory. 

Yakhot studied the equations of the structure func- 
tions by using the generating functions which had been 
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introduced by Polyakov[15, 18] . The same equations were 
obtained also by Hill by directly dealing with the velocity 
increments instead of the generating function [20]. In the 
equations, the terms arising from the convective term of 
the Navier Stokes equation are expressed in terms of 5 Pi9 
alone, but the equations again contain contributions from 
the pressure and viscous terms. Yakhot argued that in 
the equations of S2 n ,o the contribution from the viscous 
term vanishes and the pressure term can be neglected, 
and concluded C2n,o = C2n-2.2- For the equations of 
Sp.2m, the pressure term is modelled as the conditional 
average of the pressure gradient increments for given val- 
ues of U and V, (S'V±p\U, V, r), which is assumed to be 
given by the expansion in velocity increments and more 
general than Gotoh's case. He retained terms up to 
the second order, and argued a resemblance to "mean- 
field approximation" in the critical phenomena [18]. This 
model for the pressure term enabled him to close the 
equation for the probability density function (PDF) of 
V in the inverse cascading range in the two dimensional 
turbulence where the dissipation term is irrelevant. The 
PDF is found to be Gaussian and the normal scaling for 
So, q was obtained. The same closure for the pressure 
term and a model for the dissipation term yield a closed 
equation for the PDF of V in three dimensional turbu- 
lence, from which the peak value and the asymptotic tail 
of the PDF are derived, and the scaling exponent for S ,g 
is computed. 

These facts suggest that when the dissipation contribu- 
tions to the equation of the velocity structure functions 
are irrelevant the pressure- velocity correlation is a key to 
understand the anomalous scaling exponents. The aver- 
age of the pressure term conditioned upon the velocity 
increment is a kind of statistical projection of the pres- 
sure force on the dynamics of the velocity increment and 
study of its functional form would yield important knowl- 
edge of the scaling of the structure functions[18, 20]. 

Kurien and Sureenivasan studied experimentally the 
equations of the above structure functions by using the 
turbulence data measured in the atmospheric bound- 
ary layer at the Talyor microscale Reynolds number 
R\ w 10700. The width of the inertial range is about one 
decade. Since the direct measurement of the contribu- 
tions of the pressure and viscous terms are difficult, they 
considered first the equations of the structure functions 
without the pressure and viscous terms, and then exam- 
ined the balance among the terms in the equations. They 
concluded that the pressure effects are less important for 
£271,0 at low order[24]. Further, in order to assess the 
contributions of the pressure term, they applied Yakhot 's 
theory to the pressure term. With this model, the corre- 
lation between the velocity increments and the pressure 
gradient increment is expressed in terms of the veloc- 
ity structure functions alone, meaning a closure. They 
concluded that the pressure and viscous terms are small 
compared to the other terms consisting of S p . q . Also the 



scaling exponents of S'o.q, the peak values of the PDF 
and the asymptotic tail of the PDF of V were examined 
and compared to Yakhot's prediction. An agreement was 
observed for these quantities, but insufficient for the pres- 
sure contributions. Since there is small but finite degree 
of anisotropy at second and third order structure func- 
tions, and the data in the experiment suffers from noise, 
the conclusion is very suggestive but inconclusive. More- 
over the pressure and viscous terms are examined indi- 
rectly through the assumed model to the pressure term. 
Certainly we are left with lack of precise knowledge for 
the equation and further study of the equation of S p . q is 
highly needed. Especillay the direct examination of the 
pressure and viscous terms is required. 

The above story motivates us to examine the dynami- 
cal relations for S m . n (r) in terms of DNS. DNS can pro- 
vide precise data of each term in the equations so that 
the contributions of the pressure and viscous terms arc 
directly examined. When compared to the atmospheric 
boundary layer data, R\ — 460 of the DNS is low, but 
degree of isotropy is so well that about one decade of the 
inertial range width is achieved [26]. In the following, the 
analysis of the DNS data shows that the pressure term 
is important to the dynamical equations of the longitu- 
dinal correlation functions. Conditional average of the 
pressure gradient increments is also analysed and its role 
in the dynamics is discussed. 

DYNAMICAL RELATIONS 

The scalewise energy budget is well described by the 
Karman-Howarth-Kolmogorov equation 

-er = -63,0 + &v —^f- + Z force{r), (8) 

where Zf orce (r) is an energy input by an external force. 
Note that this equation does not contain the pressure 
term. When the external force has a support compact 
spectrum at low wavenumbers, both the viscous and ex- 
ternal force terms can be neglected so that we have Kol- 
mogorov's 4/5 law. 

The dynamical relations at higher order in the iner- 
tial range are derived from the Navier-Stokes equation. 
This can be efficiently done by introducing the generat- 
ing function Z = (exp(X a w a )) where w = u(x\) — u{x2) 
is the velocity difference [15, 18]. The resulting equation 
is 

dS p , q d-l + q (p-i)( g -i + d) 

~dT + r 6 ™ HjqTl) Sp - 2 ' q+2 

= -(p-l)((5V llP )UV- 2 V) 

-q((SV xp)^- 1 ^" 1 ) - D p , q , (9) 

where d is the dimensionality. D p>q denotes the contri- 
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tuitions from the viscous terms 

D p , q = „((p l)(p - 2)((V X U) 2 U?- 3 V«) 
+2q(p-l)((\7 x U) ■ (VxVpP- 2 ^- 1 ) 
+^-l)((V x ^) 2 (7 p - 1 ^- 2 )), (10) 



where X = (xi + x 2 )/2, V X U = dui(xi)/dx\ — 
dui(x 2 )/dx 2 , Vjf V = du t {x\)/dx\ - du t (x 2 ) / dx 2 , and 
ui and ut are the longitudinal and transverse components 
of the velocity vector, respectively. The external force 
term is also neglected[18, 20], and the pressure terms ap- 
pear in Eq.(9). 

In the inertial range, Yakhot, Hill and Boratav, Hill, 
and Kurien and Sreenivasan have discussed Eq. (9) for 
some set of (p, q) without pressure and viscous contri- 
butions [18-20, 24]. Following their discussion and for 
latter convenience we write a first few series of Eqs.(9) 
without the pressure and viscous terms as 



dS, 



4.0 



2. _ 6 

— 1 04. — ~ <->2,2j 

dr r r 



dSe t o 

dr 
dSs,Q 



2 _10 

~^6,0 — &4,2i 



2 14 

H 5*8.0 = ^6,2) 

dr r r 
for the longitudinal structure functions and 



dS 2 



dr 
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dr 
dS A , 2 

dr 
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'0,4, 



H — 5*2,4 — —5*0.6, 

r or 

4 4 

H — S4 = -52,4, 

r r 



(11) 
(12) 

(13) 

(14) 
(15) 
(16) 



for the mixed ones [20]. If the pressure and viscous 
terms can be neglected in the inertial range, we obtain 
the relations between the scaling exponents of the lon- 
gitudinal, mixed, and transverse structure functions as 
C2«,o = C2n-2,2, C2,2n = Co,2n+2, and so on. In the fol- 
lowing we examine the above relations quantitatively in 
terms of the DNS data. 



EXAMINATION WITH DNS DATA 

The turbulence in a periodic cubic box with size of 
2ir is simulated by using the Fourier spectral method for 
space, and Rungc Kutta Gill method for time advancing. 
The Gaussian, white-in-time random solenoidal force is 
applied to the first few shells in the wavevector space 
to maintain statistically steady state. The total number 
of grid points is N = 1024 3 and kmaxfj — 0.96, where 
77 = (v 3 /e) 1 / 4 is the Kolmogorov scale. The value 0.96 is 
slightly smaller than unity, but this does not affect the 
statistics of the inertial range. The ensemble average is 
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r/r\ 

Fig.l 4/5-law at R x = 460 [26]. 

taken over time average with the period of about four 
large eddy turnover time. Detail of the simulation is 
found in Gotoh et al.[22, 26]. 

Figure 1 shows each term of the Karman-Howarth- 
Kolmogorov equation (8) divided by er. When r de- 
creases, the forcing term decreases and — Ssfi(r) quickly 
rises to the 4/5 level and decays in the dissipation range. 
The energy budget follows Kolmogorov's 1941 theory and 
the inertial range certainly exists. 

In Figs. 2, 3 and 4, both sides of Eqs. (11), (12) 
and (13) are plotted, respectively, and the ratio of the 
difference of both sides to the left hand side (\l.h.s. — 
r\/i.s|/|?./i.s.|) is also shown in the figures. If the con- 
tribution from the pressure gradient is negligible, both 
curves must collapse at least over a part or whole of the 
inertial range. Figures show that there exist differences 
and the left hand side of Eq. (11) in Fig. 2 is bigger by 
factor 1.2 than the right hand side of Eq. (11). Dis- 
tance between the two curves is almost constant over the 
inertial range, but increases with the order n. Corre- 
sponding to this, the ratio of the difference between both 
sides to the left hand side is nearly constant over the in- 
ertial range, which suggests that the contribution from 
the pressure gradient term scales almost in the same way 

as 52 n ,0- 

When the pressure term is included in the right hand 
side of Eqs. (11), (12) and (13), the curve of the right 
hand side (dash dotted line with plus symbol) perfectly 
collapses with the curve of the left hand side, which 
means that the pressure gradient term has the same order 
of contributions as the mixed correlation term. We have 
examined the contributions from the dissipation term 
and found that they are very small compared to the other 
terms of Eq. (9) with p = 2m and q — 0, consistent with 
the theoretical prediction [18]. For the dynamical equa- 
tions of the odd order moments such as 5*271+1,0 (?"), all 
the terms would contribute equally but unfortunately we 
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have not computed these terms. 

Figures 5, 6 and 7 show the both sides of Eqs. (14), 
(15) and (16). The curves for both sides of these equa- 
tions are very close to each other for whole or part of the 
separation below the integral scale. Although it is not 
clear whether the difference between the curves remains 
to be small when the Reynolds number is increased, the 
degree of contribution of the pressure gradient to the dy- 
namical equations in these cases is weaker than the case 
of the longitudinal case. 

These observations suggest that as far as the dynamics 
of the longitudinal velocity increments is concerned, the 
pressure term is very essential, while it is secondary for 
the dynamics of the transverse velocity increments. If 
the pressure term is absent, sharp shock is formed as in 
the Burgers turbulence, but in reality the pressure sup- 
presses the formation of shock. This anti-shock action of 
the pressure would become stronger when \U\ with nega- 
tive sign becomes larger. This implies that the larger the 
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longitudinal velocity increments is, the more significant 
the contribution of the pressure term becomes (see more 
discussion in Sec. 4). On the other hand, near the shock 
region, transverse velocity increments are closely tied to 
shearing motion as results of shock killer action of the 
pressure. Both positive and negative transverse velocity 
increment occur with equal probability so that the cor- 
relation between the transverse velocity increments and 
pressure gradient would have a good chance to be small. 

When the pressure term is included in the dynamical 
equation Eq.(9), it is not clear whether or not the scaling 
exponents ^2n,o & n d C2n-2,2 are equal. We need precise 
knowledge about the correlation between the increments 
of the pressure gradient and longitudinal velocity. In the 
following we examine the pressure term from the view 
point of the conditional average. 
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CONDITIONAL AVERAGE OF THE PRESSURE 
GRADIENT INCREMENT 

The correlation function ((5p tX )U 2n ~ 2 ) can be written 
as follows 



((Sp, x )U 2 - 2 ) 



(5p, x \U,r)U 2n - 2 P(U,r)dU, (17) 



where (Sp^ x \ U, r) is the conditional average of the pres- 
sure gradient increment for a fixed value of U and P(U, r) 
is the PDF of U with separation r. If {Sp iX \U, r) is given 
as a polynomial function of U, then ((Sp. x )U 2n ~ 2 ) is 
expressed in terms of the moments of U. For the mo- 
ment we assume that (6p iX \U,r) is a quadratic function 
of U. Underlying motivation for this is that (1) the 
source term of the Poisson equation for the pressure is 
quadratic in the velocity gradient, i.e., when v! = we 
have p' = j 2 p, and (2) for large U the local Reynolds 
number is quite large so that the fluid motion would fol- 




Fig.8 — (Sp x \u, r) N for various r/rj {r m in/rj = 3 and 
rmax/rj — 1565.). As the separation r increases the 
curvature becomes smaller. 
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Fig. 9 — (5p x \u, r) N plotted against u 2 for r/r/ — 3 and 6. 

low Bernoulli's theorem locally, implying that the pres- 
sure field is proportional to the square of the velocity. 
This is sketched as follows. 

Consider a vortex tube and a cylindrical surface sur- 
rounding the tube. On this surface, the velocity vector 
is nearly along the circumferential direction and the vor- 
ticity vector is approximately parallel to the axis of the 
vortex tube. Then the vector u x oj is parpendicular to 
the surface, the Bernoulli surface[23]. On this surface, 
the Bernoulli theorem is p(x\) + u 2 (xi)/2 = Hi where 
p = 1 is assumed for simplicity and Hi is a Bernoulli 
constant on the surface. Take the gradient, we have 
~ V1P1 = U\ ■ V1K4. Then subtract from this equa- 
tion the same expression evaluated at X2 on a different 
Bernoulli surface, and introduce the new cooridnate as 
X = (xi + X2)/2 and r = xi — X2, we obtain 

-SVp(X, r) = V(X, r) ■ V x Su{X,r) 

+5u(X,r) -V r 5u(X,r), (18) 

where V(x,r) = (u(xi) + u(x2))/2 which would be a 
slowly varing function of X and r compared to Su. If 
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Fig. 10 — {Sp x \u, r) N for r/rj in the inertial range. 
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Fig. 11 Close up of — (5p x \u, r) N for r/rj in the inertial range. 

we expect \7 j^Su(X ,r) oc Su(X,r) and VrSu(X,r) oc 
5u(X,r), Eq.(18) becomes quadratic in Su. Moreover, 
since the Bernoulli constant differs on each Bernoulli sur- 
face, this fact would introduce an additional factor 5H, 
independent of the velocity difference, to Eq.(18). The 
final expression would be —5\7p = Aq + A\5u + A25u5u, 
where A/s are tensor functions of X and r. The aver- 
age over X or statistical ensemble conditional upon the 
velocity increments would yield the quadratic function of 
Su for — (SVp\Su, r). The above argument is very rough 
and needs critical examination. 

The DNS data of the conditional average at R\ = 460 
is shown in Fig. 8 for the separation from the dissipa- 
tion to the integral scale. In the plot, u and {Sp, x \u,r) N 
arc normalized by their standard deviations, respectively 
The curves are very close to the quadratic function of u. 
As the separation increases the curvature decreases, and 
the curves peel off the quadratic form at large ampli- 
tude of u. We infer that this is due to the fact that the 
number of large amplitude events becomes smaller as the 
amplitude u increases. 

To see how close the curve is to the quadratic function 
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Fig. 12 Variation of the normalized variance of the pressure 
gradient increment f-FVp e ^ 3 ^ 2 ^ 1 ^ 2 ((Sp x ) 2 ) with separation. 

of u, (5p. x \u,r) N at r/rj = 3 and 6 is plotted against u 2 
for both positive and negative u in Fig. 9. Curves are 
straight up to about u 2 w 200, strongly suggesting that 
the curves are quadratic funtions of u. The curves for the 
inertial range separation are plotted in Fig. 10. For small 
value of u, (Sp tX \u,r) N is negative and positive at large 
u. Close inspection of Fig. 11 reveals that the vertex of 
the quadratic function locates at negative u and moves 
towards the origin as the separation increases. 

It is very interesting to see how the quadratic function 
in U for {5p, x \U, r) is fixed. The assumed quadratic form 
is 

-(5p tX \U,r)=a 2 {r)U 2 + a 1 {r)U + a (r). (19) 

From the condition (Sp^ x ) — 0, we have ciq = —<i2S2.o(r) 
where we used (U) = and (U 2 ) = S2,o(r)- Also the 
condition (U8p^ x ) = by the homogeneity, isotropy and 
incompressibility yields a\{r)(U 2 ) = — ([/ 3 )a2(r)[ll, 25]. 
Combining these two relations we obtain (5p, x \u, r) N in 
non-dimensional form as 



-($P,x\u,r) N 

= -{(5p,,) 2 y 



1/2 



a-2{r) S 2 . (r) r 2 
((^x) 2 ) 1/2 r 



(Sp, x \U,r) 

[u 2 - s 3 (r)u - 1] , (20) 



where a 2 (r) = ra,2(r) > 0, and S3(r) = 
<S , 3,o( r )/ [-5*2,0 ( r )] 3 is the skewness factor. Since 
£3,0(0 = — (4/5)er in the inertial range, we obtain 



IC2.0 



(21) 



where S 2 ,o{r) = Co^^r 2 / 3 ^/^ 2 -"- 2 / 3 is used, and L is 
a macro length scale, C = %r(l/3)K w 1.32K, and K 
is the Kolmogorov constant. We expect that a 2 (r) tends 
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Fig.13 Scaled -{Sp x \u,r) N . 

to a positive universal constant in the inertial range. 
The variance ({Sp. x ) 2 ) = 2((p 2 x ) - (p, x (x + rei)p >x (x))) 
quickly approaches the value 2(p 2 x ) when r becomes 
larger than the dissipation scale as seen in Fig. 12, so that 
it is well approximated as 

<(^) 2 ) = l^v-^F^R^ (22) 

where Fy p is a nondimensional constant defined by 
((Vp) 2 ) = e 3 ^ 2 iy-^ 2 F Vp . F^ p (R x ) is R x dependent for 
low to moderate R\ and slowly tends to a constant [22, 
27]. Using these quantities we write the conditional av- 
erage as 



■(8p tX \u,r) N 
= c 2{r) 



N(r)| 



ks(r)| s 



T \ C2,0 — 1 



X 



(z) 



u + 



\ss(r)\ 



+ 1 



in the inertial range, where 

AC 



X = 



y/{2/3)Fv p (R x ) 



R 



'1/2 



(23) 



(24) 



and A is a nondimensional constant. Figure 13 shows 
the scaled plot of -c 2 (r)- 1 (5p, x |u, r) N + + l) 

against u + 



<N+{ [ 4 

for the inertial range separation. Col- 



lapse of curves for low to moderate amplitude is excellent. 

We have computed the curvature of —(6p tX \u) N at the 
origin and plotted in Fig. 14. Since c 2 (r) is assumed to 
obey the power law c 2 (r) = x(rjV)~ a in the inertial 
range, the curve is multiplied by (r / L) 1 ^'' 2 - with (2,0 — 
0.701 (see Fig. 14) [26]. The nondimensional constant \ 
and the exponent a were found to be 0.027 and 0.267, 
respectively, by the least square fit in the range of 0.09 < 
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Fig. 14 Compensated plot of the coefficients C2(r). 
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Fig. 15 The role of the pressure term on the dynamics of the 
longitudinal velocity increment. The convection stretches 
the fluid portion with 8U > (£ = dll/dx decreases) and 
squeeze the fluid portion with SU < (|£| increases). The 
modeled pressure gradient gives positive contribution C2U 2 

to dU/dt, that is, £(> 0) of the stretched portion is 
increased and |£|(£ < 0) of the squeezed one is decreased. 
This means that the stretched and squeezed portions are 
pushed back left by the pressure gradient. 



r/L < 0.3. The latter value 0.267 is slightly smaller than 
the value 1 - (2,0 = 1 - 0.701 = 0.299. 

From the above observation, we conclude that the con- 
ditional average — (Sp iX \u, r) N is well approximated by 
the quadratic function of u with the curvature which de- 
creases with the separation. The quadratic function of u 
for — (Sp tX \u, r) N is consistent with the anti-shock action 
of the pressure. Locally the equation for U may be de- 
scribed by the Burgers type equation with the pressure 
effect: 



DU 
~Dt 



dU 
~dt' 



or 



dS P d 2 U 
~d~X + "Ifr 2 



■K(X,r,t), (25) 
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where X is a x-coordinate moving with the average veloc- 
ity at two points x\ and x 2 , and K denotes the remain- 
ing terms arising from the Naveir-Stokes equation[14]. If 
the pressure gradient term is modeled by the quadratic 
function of U as inferred from the conditional average 

— (dp,x\U,r), all the terms except K term are expressed 
in terms of U . This replacement is thought to be a kind 
of statistical projection of the pressure gradient onto U- 
space. The K may be modeled as random force, but it is 
irrelevant in the following discussion. 

To make discussion simpler, we take the leading term of 

— (Sp, x \U,r) « a 2 (r)U 2 ,a 2 > 0. In the large negative U 
region (squeezing domain), the convective term produces 
a shock as in Fig. 15. The pressure gradient term makes a 
positive contribution a 2 U 2 to dU/dt, which means that 
large (negative) amplitude U(t + dt) is decreased. On the 
other hand, for large positive U regions, those regions are 
stretched by the convective action so that the positive 
amplitude decreases[12, 13]. When the pressure gradient 
term is added, again this gives positive contributions to 
dU/dt, thus U(t + dt) increases. Therefore the pressure 
gradient term acts in such a way that it suppresses the 
formation of singularity such as shocks by the convective 
motion, and this can be understood as the action which 
the curve U(t + dt) is pushed back toward left by the 
pressure gradient in Fig. 15. This implies that the pres- 
sure term is an indispensable ingredient for the dynamics 
of the longitudinal velocity increment. 



IMPLICATION OF THE PRESSURE GRADIENT 
TERM 

When the closure (19) is included in Eq.(9) withp = 2n 
and q = 0, the pressure- velocity correlation term becomes 



2n-l,0 



(r) 



S2,o(r)S2„-2,o(r)Y (26) 



Suppose that when the Ryenolds number and n are very 
large, S 2n fi{r) m Eq.(26) is dominant in the inertial 
range. This requires ( 2nfi < C2n-i,o + 1 - (2,0- On the 
other hand, the realizability condition requires that Cp,o 
is a nondecreasing function of p, meaning £271-1,0 < C, 2n ,o- 
Both constraints suggests that < ( 2n ,o — C2n-i,o < 
1 - C 2 ,o w 0.3. In K41, C 2 „,o - C2n-i,o = 1/3- From 
the DNS data, the above requirement is satisfied about 
at the order greater than four. 

Dominance of S 2n ,o over other two terms of Eq.(26) 
means that Eq.(9) becomes 

dS 2n , 2 _ 2(2n-l) 

; H-O2n,0 — J2n-2.2 + \^n— 1) O 2 n,0, 

dr r r r 

(27) 



asymptotically in the inertial range. Suppose that 
d 2 (r) cx r 7 in the inertial range. When 7 < 0, since 
the two terms of the right hand side of Eq.(27) are posi- 
tive and the last term becomes dominant, the both side 
can not balance. This means that the case of 7 < can 
not occur. When 7 = 0, both sides of Eq.(27) balances, 
thus we have C2«,o = C 2 n- 2 .2- If 7 > 0, the last term of 
the right hand side of Eq.(27) becomes less important so 
that the left hand side balances with the first term of the 
right hand side and C2n,o = C2«-2,2- The conclusion of 
the above arguments is C2n,o = C2n-2,2 for large n, but 
this equality is totally dependent on the assumption of 
the quadratic function in U for the conditional average 
(see also the discussion by Hill [20]). 

Although the DNS data shows that d 2 (r) = 
xC _1 (2Fvp/3) 1 / 2 is approximately constant, there is 
small but finite difference between the exponent a = 
0.267 and 1 — (2,0 = 0.299. If we regard these two values 
correct, we have 7 = 0.299 - 0.267 = 0.032 > 0. How- 
ever, we think that our data is not accurate enough for 
definite judgement, and the preceding argument about 
the role of the pressure term suggests that 7 = is very 
plausible. 

It is important and interesting to seek physical ex- 
planation for possible r-dependence of a 2 (r), if it is the 
case. To the authors' best knowledge, no theory has been 
developed to compute a 2 (r) including a nondimensional 
constant. We present one possible argument in the fol- 
lowing. 

The pressure is given by the spatial integral of the 
source term which is quadratic in the velocity gradients, 
so that the pressure gradient is given by 



Vp=- I V x G(x-y)B(y)dy, 
Jv 



G(x - y) 



Air \x — y\ ' 



(28) 
(29) 



for the periodic boundary condition in the present study, 
where B = Vm : Vm. The integral is roughly evalu- 
ated as the product of B, VG and an effective volume 
V e f feet [28] . The derivative of the Poisson kernel VG does 
not yield non-integer power, and the source term B would 
contribute to S 2n _o. The effective volume is the one over 
which the source term has a dominant contribution to 
the integral, and would change its shape and size de- 
pending on the amplitude and its spatial coherency of B. 
This effective volume is characterized by a geometrical 
non-integer dimension. Therefore it is plausible that the 
effective volume V e f f ec t may induce a non-integer power 
law dependence of a 2 (r) on the scale r. 

SUMMARY 

We have examined the dynamical relations for the var- 
ious kinds of the velocity structure functions in terms of 
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the DNS data. It was argued that the pressure term is 
a key to understand the scaling of the velocity structure 
functions. The contribution of the pressure gradient is 
important to the dynamics of the longitudinal structure 
functions, while in the case of the transverse or mixed 
ones the pressure gradient may or may not be important. 

To assess the statistical role of the pressure gradient 
term in the dynamics of the longitudinal velocity incre- 
ment, we have computed the conditional average of the 
pressure gradient increment with a given value of U and 
r. It was found that the quadratic function in U is a good 
approximation to the conditional average. By applying 
two constraints to the conditional average and using Kol- 
mogorov's 4/5 law, we have theoretically determined two 
coefficients of the quadratic function but there remains 
one undetermined coefficients. 

The role of the pressure gradient was argued in terms 
of the equation for U which is similar to the Burgers' 
equation but with the term which models the pressure 
gradient. The pressure gradient acts to resist to stretch- 
ing and squeezing of the fluid element, leading to shock 
smearing, which in turn means that the scaling exponent 
of 5271,0 is n °t constant in n, unlike the Burgers case, and 
is a slowly increasing function of the order. This means 
that the pressure is an intermittency killer (Kraichnan 
1991) [29]. 

The quadratic function form for the conditional aver- 
age of the pressure gradient increment leads to the closure 
of the moment equation for S2 n ,o- If higher order poly- 
nomial form is used, the closure fails. However, if a thin 
vortex tube is the most singular object under the Navier- 
Stokes dynamics, the argument using the Bernoulli the- 
orem nearby the tube described in Sec. 4 suggests that 
there are no higher order terms than the second order. 
The bottom line is C2n.o — C2n-2,2- It was argued that 
the effective volume is an important factor for the de- 
pendence of &2(r) on r. Further theoretical study is nec- 
essary. The conditional average of the pressure gradient 
also leads to the closure of the Liouville equation for the 
PDF of U. This analysis is now underway and will be 
reported elsewhere. 
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